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Abstract. We consider a cable described by a discrete, space-homogeneous, 
quasi one-dimensional Schrodinger operator Hq. We study the scattering by a 
finite disordered piece (the scatterer) inserted inside this cable. For energies E 
where Hq has only elliptic channels we use the Lippmann-Schwinger equations 
to show that the scattering matrix and the transfer matrix, written in an 
appropriate basis, are related by a certain polar decomposition. For energies 
E where Hq has hyperbolic channels we show that the scattering matrix is 
related to a reduced transfer matrix and both are of smaller dimension than 
the transfer matrix. Moreover, in this case the scattering matrix is determined 
from a limit of larger dimensional scattering matrices, as follows: We take a 
piece of the cable of length m, followed by the scatterer and another piece of the 
cable of length m, consider the scattering matrix of these three joined pieces 
inserted inside an ideal lead at energy E (ideal means only elliptic channels), 
and take the limit m — > oo. 



1. Introduction 

We consider discrete quasi one-dimensional Schrodinger operators on strips of 
width N of the form 

(1.1) (i^*)„ - -*„+i - + 

where * = (*„)„ez e ^^(Z, C^) ^ e^{Z) is an sequence of vectors in 

and Vn E Hcr(Af) is a bounded sequence of Hermitian N x N matrices. Such an 
operator is a so called tight binding model for a cable with N channels. The terms 
— \E'„+i — Vl/n-i correspond to the horizontal Laplacian and describe the 'hopping' 
of an electron from state to state along the wire. The matrix potentials Vn describe 
the hopping or interaction between the different channels and may also include some 
potential. A particular case of interest are models where the Vn are perturbations 
of a fixed matrix W. li Vn = W for all n then one finds Bloch waves and the 
operator describes a pure, space homogeneous cable with a pure crystal structure. 
The perturbations then model impurities in the cable. For instance, randomly 
doped semiconductors are supposed to be modeled by random potentials Vn- For 
instance, the case where A'^ = 1 and the Vn are independently identically distributed 
corresponds to the one-dimensional Anderson model as proposed by Anderson J]. 
Choosing Vi to be distributed according to the Gaussian unitary ensemble (GUE) 
or the Gaussian orthogonal ensemble (GOE) for N > 1 corresponds to a Wegner 
iV-orbital model. WegnerjTT] studied the N oo limit of such models. 

The focus in mathematical physics often lies in the spectral theory on the infinite 
strip. From a solid state physics point of view the electronic properties of finite 
pieces are quite of interest. The general idea is that absolutely continuous spec- 
trum corresponds to a conductor even for infinitely long pieces, whereas Anderson 
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localization corresponds to an isolator when the length of the piece is much larger 
than the localization length. 

In the theory of electronic conduction as developed by Landauer [5J [51 , Imry [7] 
and Biittiker [H [5] a scattering approach is used. The idea is that the electronic 
properties of such a finite piece, from now on called the scatterer, is in principle 
given by considering the scattering of this piece inserted inside an ideal lead. By 
ideal lead one means a pure cable with only elliptic (propagative) channels. The 
mathematical definition will be given in the next section. The scattering matrix 
for this scattering problem describes reflection and transmission of incoming Bloch 
waves to outgoing Bloch waves on the right and left of the scatterer. Related in a 
twisted way to this scattering matrix is the S -transfer matrix giving the transfer 
from waves from the left to the right of the scatterer. We call it S'-transfer matrix 
as it is obtained from the scattering matrix and we will use the terminology transfer 
matrix for a different object. 

From the scattering matrix or the S'-transfer matrix one can calculate certain 
quantities such as the Landauer conductance or shot noise. For more information 
on these connections I recommend the review by Beenakker [3] . The scattering and 
the S'-transfer matrices depend on the specific choice of an ideal lead as well as 
on the choice of a basis for its Bloch waves, but important quantities such as the 
Landauer conductance do not. 

The advantage of the S-transfer matrix compared to the scattering matrix is the 
so called multiplicity property. The physics intuition is the following. Suppose one 
puts two scatterers together which are described by S-transfer matrices Ti and 
Then the first transfer matrix Ti connects the amplitudes and phase information 
of waves on the left of scatterer 1 to the right of scatterer 1 which is the left of 
scatterer 2. Now, T2 connects these amplitudes and phases to the ones on the right 
of scatterer 2. Therefore, the product T2T1 connects the amplitudes and phases 
on the left of the two scatterers to the right of the two scatterers. Thus, T2T1 
corresponds to the S-transfer matrix of both pieces put together. 



In the mathematical analysis of operators as given by (2.1 ) one defines the trans- 
fer matrix from the stationary Schrodinger equation (cf. (2.2 1 and (2.4)). These 
transfer matrices satisfy the multiplicity property which can be seen easily. For an 
ideal lead the transfer matrix is conjugated to a unitary matrix. If one diagonalizes 
it then it looks like the S-transfer matrix of Bloch waves. In fact, it seems to be 
quite known that using the same basis change of a disordered piece corresponds to 
the S-transfer matrix of this piece w.r.t. the same ideal lead. For instance, this is 
mentioned in Ref. [2] and it will be confirmed in this article. 

An important development in the electronic conduction theory is the so called 
DMPK[6l|T2] theory and DMPK equation. This is a stochastic differential equation 
(SDK) describing the conductance of a disordered wire with respect to its length 
in a macroscopic setup. Bachman and de Roeckp^ analyzed the connection of the 
microscopical Anderson model on a strip to DMPK theory. If the unperturbed 
operator describes an ideal lead, then they found an SDE describing the evolution 
of the transfer matrices in an appropriate scaling limit. This can not be obtained if 
the unperturbed operator is a pure cable with elliptic (propagative) and hyperbolic 
(non-propagative) channels. I believe that in this case one should consider the 
S-transfer matrix coming from scattering a disordered piece with respect to the 
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unperturbed operator. In these cases the scattering matrices and the S'-transfer 
matrices are of lower dimensions than the transfer matrices. Also, the multiplicity 
property for the S'-transfer matrices is no longer valid, but it still holds for the 
transfer matrices. The purpose of this paper is to analyze the relations between 
these matrices in this case (cf. Theorem 2.1 1. 

From a physics point of view, the scattering matrix of a finite disordered piece 
with respect to a cable with hyperbolic channels does not only contain information 
about the scatterer but also about the cable. This is in principle also true if one 
has an ideal lead, but since an ideal lead has only propagative channels, it does not 
affect important quantities such as the Landauer conductance. However, hyperbolic 
channels do have an effect. Therefore, it should be treated as a scatterer itself. 
By physics intuition, the situation of having the finite scatterer inserted inside an 
infinite cable should be described by the following limit: We take a piece of the 
cable of length m, followed from the finite scatterer and another piece of length m 
of the cable., Then we obtain the scattering matrix for these three blocks together 
inserted in an ideal lead and take the limit m ^ oo (cf. Figure [l] on page [7|. We 
will prove that this limit gives indeed the scattering matrix of the scatterer with 
respect to the pure cable with hyperbolic channels, cf. Theorem |2.3[ 



Acknowledgment: I am thankful to H. Schulz-Baldes and A. Klein for many 
suggestions. 

2. Statement of Results 
As described above, let H be an operator on ^^(Z,C^) = £'^(Z) (8)C^ defined by 

(2.1) (F*)„ = - + K^-™ 

where Vn G Her(A'') is a sequence of Hermitian N x N matrices. H describes a 
cable with A'^ channels. Associated with such an operator are the transfer matrices 
■ They arise from the stationary Schrodinger equation H'i> = E'i', which gives 

Note that T^f is in the conjugate symplectic group Sp(2A^) defined by 

(2.3) Sp{2N) = {Te Mat(2Ar, C) : T*JnT = Jn} where Jn = 

The individual blocks are all of size NxN. This group is different from the complex 
symplectic group Sp(2Af,C) = {T : T'^ JnT = Jn}- 

The transfer matrix of the block of length L — I from ^ to L — 1, where / < i, is 
given by the product 

(2.4) 71^. ^rf_irf_,...7I^ which gives T^l = 

if ijvj/ = E'^. This product only depends on E and the sequence V;, . . . ,Vl-i- 
Hence, for a fixed energy E, each such sequence gives rise to a certain transfer 
matrix. Moreover, the transfer matrix for two consecutive blocks (sequences), is 
just the product of the transfer matrices for each block, e.g. 7^^^ = 7^^ 7^^ for 
< I < L. We referred to this as the multiplicity property in the introduction 
above. 
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We want to insert a finite block of lengtli L witliin a space-homogeneous cable 
and consider it as a scatterer within the cable. The scatterer will be described 
by the sequence Vb, . . . , Vl-i and the transfer matrix 7^^^ which connects to 

* solution of imi = E"^. The space-homogeneous cable will be described 
by the operator 

(2.5) (i/o*)n = - *„-i + W-^n , * = (*„)„ez e £2(z,C^) . 

The difference to H is that the Hermitian matrix W is always the same and Hq is 
invariant by translations on Z. Therefore we call Hq space-homogeneous. 

Inserting the finite scatterer is described by changing the operator Hq on a finite 
piece. Therefore, let {Vn)n satisfy 

(2.6) Vn^W for n < and n>L, 



then the operator H as defined in (2.1 1 describes the cable with the inserted scat- 
terer given by the sequence Vq,. . . , Vl-i- The scattering of this piece is described 
by the unitary scattering operator S = of H with respect to Hq, where 

r2± = s — limt_j.zpoo e'*^e~''^''. Since S commutes with Hq, it can be represented 
by scattering matrices on the energy shells for (almost) each energy E in the spec- 
trum of Hq . 

The Hermitian matrix W describes the transverse modes in the cable. Let tpa S 
C^, a = 1, . . . , be an orthonormal basis of eigenvectors of W with corresponding 
real eigenvalues Aq. The spectrum of Hq is purely absolutely continuous and given 
by the union of N bands, Ua=i[-2 -I- A^,, 2 -I- A^]. Given an energy E, ipa is called 
an elliptic channel if |Aq — i?| < 2, a parabolic channel if |Aq. — i?| =2, and a 
hyperbolic channel if |Aq, — 2| > 2. If there is a parabolic channel then E is called 
a band-edge. The number of elliptic channels at E will be denoted by s{E), the 
band-edges are exactly the discontinuities of s{E). If E is not a band-edge, then 
the multiplicity of the spectrum of Hq at E is given by 2s{E) which exactly equals 
the number of eigenvalues of modulus 1 (counted with multiplicity) of the transfer 
matrix 

(2.7) -^1 

Since the multiplicity is 2s{E), the scattering matrix describing the scattering op- 
erator on the energy shell has to be a 2s{E) x 2s{E) matrix. The corresponding 
extended states of Hq can be split into s{E) right-moving and s{E) left-moving 
waves at energy E. In the sequel we will often use s instead of s{E). 

The terminology transfer matrix also appears in the scattering theory of elec- 
tronic conduction as developed by Landauer [51 [S] , Imry [7] and Biittiker [?J [S] . A 
short overview is given within a review by Beenakker [3]. We will call this transfer 
matrix the ^-transfer matrix in order to distinguish it from the transfer matrix 
as defined above. The ^-transfer matrix connects waves on the left to waves on 
the right of the finite scatterer, whereas the scattering matrix, let us call it , 
relates incoming and outgoing waves. For the scattering matrix we choose the 
following convention. Writing = ( ;^ ) , the s x s matrices T, T' correspond to 
transmission of waves from left to right, resp. right to left, and R and R' correspond 
to reflection of waves on the left, resp. right of the scatterer. Then, one has the 
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following relations, 



— I 5 



where a+ , 6+ G C are vectors describing the amplitudes of right-moving waves on 
the left, resp. right side of the scatterer, and a~,6~ G describe the amplitudes 
of left-moving waves on the left resp. right side. 

As the scattering operator is unitary, the scattering matrices are unitary as well. 



i.e. iS S U(2s). From the relation (2.8) one then finds that T is in the pseudo- 



unitary or Lorentz group U(s, s) of signature (s, s), defined by 

(2.9) U(s,s) = {f e Mat(2s,C) : f*asf = ^,,} where Gs = \ 

The blocks in Qg are all of size s x s making it a 2s x 2s matrix. The conjugate 
symplectic group Sp(2s) and the Lorentz group U(s, s) are related by the Cayley 
matrix, 

(2.10) C,Sp(2s)C: = U(s,s) , where = ^ Q G U(2s) . 

As described in Ref. [TTl [T^l |3] and in Appendix [Xj and are related by 
the polar decompositions 



(2.11) = 

(2.12) = 



Ur.+ \ f VQ VQ^\ fUi,+ 



ui_ \ /-v/i-Q-1 a (Ui,+ 



Here Q is a real, diagonal matrix satisfying Q > 1 and Ui^±, Ur,± G U(s) are unitary 
matrices mixing the channels on the left and the right. As shown in Appendix \K\ 
for any pseudo-unitary matrix one finds a unitary matrix satisfying ( |2.8[ ) 
by these polar decompositions. However, given one may not always find T^, as 
the matrix < 1 occurring in the polar decomposition of is not necessarily 
invert ible. 

We will show that the iS-transfer matrix and the transfer matrix are related. In 
fact for energies where Hq has only elliptic channels, they are simply related by a 
conjugation. This is a well known fact and appears e.g. as a Lemma in Ref. [2]. 
Using the Lippmann-Schwinger equation, it will be confirmed once more. The 
new investigation in this paper is the relation if the background operator Hq has 
hyperbolic channels. Then the S'-transfer matrix is of smaller size than the transfer 
matrix and the relation between them is more complicated. More precisely, we 
obtain the following. 

Theorem 2.1. (i) There is a unitary operator U : f{Z,C^) J® C'^'^'-^^dE, only 
depending on Hq, and, except for finitely many energies in the spectrum of Hq, there 
exists a scattering matrix S lJ{2s{E)), such that the spectral decompositions of 
Hq and the scattering operator S are given by 



(2.13) Hq = U* 



U , S =U* WVl 



5^ dE 



V,U. 
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with 
(2.14) 



Vr 















AE 



Here, ks is a real, s{E) x s{E) diagonal matrix and its entries are the wave numbers 
for the extended states of Hq at energy E. 

(ii) Let Hq have only elliptic channels at E. Then as in (i) and the S -transfer 
matrix defined by (2.8) both exist. Moreover, there exists M G Sp(2iV) [defined 
in (4.3)) only depending on E and Hq, such that 



(2.15) 



N 



In particular, S and Tqj^ are related by the basis change (2.151 and the polar 



decompositions (2.111, (2.12) 



s > 
Moreover, 



(iii) For all but finitely many energies E where Hq has s > elliptic and N 
hyperbolic channels there exist as in (i) and defined by (2.1 
there are matrices M G Sp(2A^) [defined in (4.3)) only depending on E and Hq, 
such that 
(2.16) 



= CsOM-%^^M 



1-0* { OM-%^^MO* ) OM-'Tfl^^M 



o*c: 



whe 



O = (o 1 o) e Mat(2s X 2iV) and 6 = [0100) E Ma.t[N - s x 2N) . 

The rows are divided in 4 blocks of sizes s, N — s, s, N — s and 1 always denotes 
a unit square matrix. As the conjugation with O reduces the dimension, we call 
ClT^Cg e Sp(2s) and also G U(s, s) itself a 'reduced' transfer matrix. 

Remarks. 1. The unitary operator U is chosen such that the first s[E) entries of 
[U'^)[E) correspond to right moving waves and the other ones to left moving waves. 
The off diagonal block structure appearing in the direct integral in the definition of 
W interchanges right and left moving waves and is necessary in the used convention 
for the scattering matrix S^, as the diagonal blocks correspond to reflection, not 
transmission. 

2. The expressions g"**^^^ appearing in the definition of Vl correspond to different 
phase normalizations for waves on the right and the left of the scatterer. This way, 
if if = Hq and hence 5 = 1, then one has 



5^ = 











AkEL 











Therefore, the 5'-transfer matrix for a piece of length L of the cable Hq gives 
precisely the phase evolution of the waves in that piece. This is a reasonable 
convention for the 5'-transfer matrix. 

3. The finitely many energies, where does not exist, consist of the band edges 



(discontinuities of s[E)) and the energies where O* M~^Tq^]^MO as in (2.16) is not 
invertible. The latter is the case if either E is an eigenvalue of H or if some waves 
of Hq with energy E are totally refiected. We show that these cases happen only 
at finitely many energies E. 

In the theory of electronic conduction developed in Ref. [HO [71 [HI IS], the scatterer 
is connected to so called ideal leads. In this case the ^-transfer matrices coming from 
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scattering theory are supposed to have the multiphcity property, i.e. the iS-transfer 
niatrix for two consecutive blocks is just the product of the ones for each individual 
block. In fact, for energies E where Hq has the maximal possible multiplicity, i.e. 
s{E) = N, the S'-transfer matrix is related to the transfer matrix 7q^ by a 



simple basis change as given in (2.15). Since the multiplicity property mentioned 
above is obviously true for the transfer matrix, it follows for the S'-transfer matrix 
in this case. 

However, if Hq has hyperbolic channels, s{E) < N, then the S'-transfer matrix 
as in (2.161 does not have this property anymore. For that reason, we make the 



following definition. 

Definition 2.2. Let Hj be an operator on ^^(Z,C^) given by 

(2.17) {Hi^j)n = -^n+l - + Wjl^n , 

where Wj G Her(A^) is a Hermitian N x N matrix. Hj is called ideal at an energy 
E iff all channels are elliptic for that energy. Equivalently, this means that the 
multiplicity of the spectrum of Hj is equal to 2N in a neighborhood of E. 

If Hq is not ideal at E then one might consider it as a scatterer with respect to 
an ideal lead at the energy E. In particular, from physics intuition, the scattering 
matrix of the finite block with respect to the non-ideal lead Hq should be described 
by the following limit: Take a piece of the cable described by Hq of length to, fol- 
lowed from the finite scatterer described by the sequence Vq, . . . , and another 
piece of length to of the cable described by Hq, connect them to an ideal lead on 
the right and the left, calculate the scattering matrix and take the limit m — > cx) 
(cf. Figure [l]). 
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Figure 1 . Scatterer and pieces of cable inserted inside ideal lead. 



Therefore, let E be an energy where Hq has s < N elliptic channels and where 
the S-transfcr matrix and the scattering matrix exist. We construct an ideal 
operator Hj at energy E by defining an appropriate hermitian matrix Wj. The 
spectral decomposition of W is given by W — '^^^ ^afa'fia- Assume that the ipa 
for a < s are precisely the elliptic channels, then we define 



(2.18) 



Wi 



a=l 



N 

E 

a— s+l 



and let Hj be given by (2.17). Furthermore, let the operators H{in) be defined by 

(2.19) (i?(m)^')„ = -*„+i - ^'„_i + K(m)*„ 

with Vn{m) = V„ for < n < i— 1, Vn{m) = W for —to < n < and L < n < L+m 
and Vn{m) = Wi for n < —to and n > L + m (cf. Figure [T]). The scattering of the 
operators H{m) with respect to Hj at energy E is described by scattering matrices 
Sf{m) as in Theorem 2.1 (where H gets replaced by Him), and Hq is replaced by 



Hi). Sj{m) exists as Hj is ideal at E. 
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Theorem 2.3. Let the 2s x 2s scattering matrix at E for the scattering operator 
of H with respect to Hq be given by — ( ) , written in s x s blocks. There is 
a diagonal, unitary matrix D, and there is a real diagonal N — s x N ~- s matrix 0, 
such that for the scattering matrix Sf {m) describing the scattering of H{m) with 
respect to Hi one has 



(2.20) 



lim D'^Sf{m)D"' = 



R T' 

-e'" 

T _R' 

V e'" 



The whole matrix has size 2N x 2N and is divided in blocks of sizes s, N — s, s, N - 



In particular, using the 2s x 2N matrix O as in Theorem 2.1 (Hi), we obtain 



(2.21) 



= O 



lim D"'Sf{m)D' 



O* 



Remarks. 1. The unitary, diagonal matrices Z?™ are just phase normalizations 
counteracting the phase evolution on the pieces of length m from the cable de- 
scribed by Hq. The e*^ terms correspond to a total reflection in the hyperbolic 
channels of Hq in the limit with some specific phase change. 

2. The transfer matrix for the inserted piece in the ideal lead Hj, described 
by H{m), is given by {T°'^r'%%{T°''^)'^ which is related to Sf{m) by Theo- 
rem 2.1 (ii). In this sense, the reduced transfer matrix which is related to can 
be interpreted as some sort of limit of (T°'^) 7o^(T°'^) for m — > oo, combined 
with a projection on the elliptic channels. 

One of the interesting byproducts of this work is the reduced transfer matrix 



and its relation to the transfer matrix as given by (2.16). For the reduced transfer 



matrix, the hyperbolic channels get eliminated in a specific way. Let me briefly 
explain with some conjectures why I believe this object is of further interest. 

Assume the matrix potentials Vn are random perturbations of W , i.e. Vn = 
W + \Wn where A is small and the Wn are i.i.d. random Hermitian matrices 
with mean zero. Then the transfer and scattering matrices are random. If Hq has 
only elliptic channels, Bachmann and de Roeck [5] as well as Valko and Virag [TS] 
obtained a stochastic differential equation (SDE) for the evolution of the transfer 
matrix in the limit A — ^ 0, L = cA~^ — >■ oo. In the presence of hyperbolic 
channels, such a result can not be obtained. The main motivation for Bachmann and 
de Roeck [2] was to investigate the relation of such models to DMPK[SJ [T^] theory 
which studies transport in disordered wires using scattering matrices. Therefore, 
the reduced transfer matrix may be of interest. 

Conjecture 1: The evolution of the random reduced transfer matrix can he de- 
scribed by an SDE in the appropriate scaling limit A — > 0, L = cA^^ — > oo. 

Related to Ref. |2] and jl6] is the perturbative calculation of the invariant 
measure of the random action of the transfer matrices on the flag manifold in the 
limit A — > 0. This action is studied to obtain the Lyapunov exponents. For energies 
where Hq has only elliptic channels, Sadel and Schulz-Baldes 15] showed under 
generic conditions on the randomness, that the weak-* limit of the invariant measure 
exists and has a smooth density with respect to a canonical Haar measure. This 
weak-* limit distribution could also be obtained from the limit SDE. In the presence 
of hyperbolic channels, such a limit distribution should exist and be supported on 
a certain stable submanifold determined by the hyperbolic channels, as explained 
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by Romer and Schulz-Baldes [M] who did some numerical calculations. The stable 
submanifold is isomorphic to a flag manifold on which the reduced transfer matrices 
act. Therefore, I have the following conjecture. 

Conjecture 2: The perturhative invariant measure on the stable submanifold is 
related to a limit SDE as in Conjecture 1. 

Let me give a short outline. In Section [3] we will consider the spectral decom- 
position of Ho and scattering states of H. In Section [3] we obtain some normal 
forms of the transfer matrix after certain basis changes. Guided by physics 
intuition we define a reduced transfer matrix in Section [5] In Section |6] we use 
the Lippmann-Schwinger equations to obtain Theorem |2.1[ Finally, we show The- 
orem [23] in Section 

3. Channels and scattering states 

We will use Dirac notations, hence expressions like denote vectors ^ € (C^)^ 
(not necessarily in and for n G Z we denote the vector '^n G by (?^|^'). Let 
|n,0 for n G Z,; e {1, . . . ,iV} denote the (.'^{'L,C^) vector defined by {m\n,l) = 
Sm.nei, where ei is the l-th canonical basis vector in C^. 

As above, let ipa G , a = 1, . . . , TV be an orthonormal basis of eigenvectors 
of the Hermitian matrix W and denote the corresponding eigenvalue by A^, i.e. 
Wifa — Xafa- Furthermore, define |5'Q,fc) by 

(3.1) (nl^'lk) = ^^e^^-, 
then one finds 

(3.2) Ho\^l,k) = {-2cos{k) + X^)\^>l,k) . 

These pseudo-eigenvectors form a partition of unity in the sense that 

N ^ 

(m,j|^'°,A:)(«'°,fc|n,Odfc = S^^^S^j = {m,j\n,l) 



1 



^ 27r 



Therefore, in a weak operator topology induced by the functionals B i— >■ (m, j V) 
(as I am not testing with all vectors this topology is actually weaker than the 
usual weak operator topology) one can write 

(3.3) 1 = E^/ l*°,fc)(*a,fc|dfc- 

Extending the Fourier transform i^((— tt, tt), ^) — > £^(Z) to distributions, one 
formally obtains from the inverse transform 

(3.4) ^(.pO,fc'|M.o,fc) =Y.^^*^^e^^^-^'> = 5^j,Sik^k'). 

In this sense, the pseudo-eigenvectors form an orthogonal system. 

Recall that we called an eigenvector (fa of W an elliptic channel for the energy 
E iS \E — \a\ < 2. In that case there exists ka G (0,7r) such that 

(3.5) £; = -2cos(fc„) 4- A„ . 

The terminology elliptic comes from the fact, that this corresponds to eigenvalues 
g±jfea ^Yie transfer matrix T"*''^ and is therefore related to a rotation. Now 
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consider fc^ as a function ka{E), where the interval on which this function is defined 
depends on a. To change the normahzation of the pseudo-eigenvectors with respect 
to energy, define for the eUiptic channels 
(3.6) \K,E,±) = (4^sin(fc„))-i/2|^o^^^^^ 



which by (3.2 1 and (3.5) are pseudo-eigenvectors of Hq with energy E. A change 
of variables in p.3[ ) shows 

(3.7) 1 = E / {\K,E,+}{^lE,+\ + \^lE,-}{^lE,~\)dE, 
and the spectral decomposition of Hq is given by 

(3.8) i?o = E / E {\^l,E,+){^<i,E,+\ + \^lE,-}{^lE,-\)dE. 

Furthermore, we say that ipa is an hyperbolic channel iff |i? — A^j > 2 and a 
parabolic channel iff — Aq,| = 2. The parabolic channels correspond to band edges 
and there are at most 2N of them. Now let E be some energy in the spectrum of 
Hq without any parabolic channel. Then there is at least one elliptic channel for E. 
Let us reorder the channels such that ipi, . . . ,ips are elliptic and tps+ii ■ ■ ■ , are 
hyperbolic channels. Furthermore, for the hyperbolic channels a > s define 7q, > 
and Ua G {—1, 1} such that 

(3.9) = -2MaCOsh(7„) A" , (a > s) . 
Then the 2N vectors 

i!^" ),l<a<s, and ( ^" ] s <a<N 

are eigenvectors of 7"°'^ and form a basis of C^^. Any formal eigenvector of 
Hq satisfying Hq'^'^ — E"^^ is uniquely defined by Xl/g and ^P]* and hence a linear 
combination of the 2N formal eigenvectors \^^^,E,+), \^^^,E^—) [a < s) and 
\i'l,E,+), 1*°,^;,-) («>s) given by 

(3.10) (n|*°,i;,e) =(p,[27rsinh(7„)]-5uie-+<e®^"", a>s,ee {+,-}. 

Here and below we use as variable symbol for + or — . In this sense, 5e,+ = 1 
for Q — + and (5e,+ — for = — . For a number c we define 0c = ^e,+c — 6q^-C. 
The factor in (3.10) seems strange but it leads to nice relations in the next section. 

Thus, for a formal eigenvector |^'°,£') of H^ there are coefficients C^jC^ for 
a < s and , for a > s such that 

{n\^",E) = J2 + E ^e|#°,£;,0) 

Q<s,ee{+,-} Q>s,ee{+,-} 

Now let I^I^, E) be some formal eigenvector of H with eigenvalue E. Then for n < 
and n > L — 1 it looks like a formal eigenvector of Hq. Therefore, there are constants 
a+,a~,6+,6~ for a < s and aJ,a~,6+,6~ for a> s associated to by 

(3.11) {n\^,E)= E a2{n\^l,E,Q) + E &2{n\^l,E,e) , 

a<s a>s 

ee{+-} e6{+-} 



RELATIONS BETWEEN TRANSFER AND SCATTERING MATRICES 



11 



for n < and 
(3.12) 

a<s a>s 

ee{+,-} eG{+-} 

for n > i — 1. 1^1', £■) is an eigenvector of H iff there are only exponential decay- 
ing parts for the limits n — > ±c», which means that a'^ = = = = 
0, = 6+ = where a"*", a"*" denote the vectors (a+)i<Q<s, ia'^)s<a<N and 
a^, a^, 5+, b^ ,b^, b~ are correspondingly defined, j^*, E) is called a scattering state, 
extended state, or pseudo-eigenvector of H iff it is not an eigenvector and has no 
exponential growing parts, neither at -|-oo nor at —oo, which means — b'^ — 0. 
(These are the states that can be used to create a sequence \'^n) of normalized 
£^ vectors by cut offs, such that \\{H — E)\'^n)\\ for n oo. Hence by the 
Weyl criterion, E is in the spectrum of if a scattering state exists.) Thus, a 
pseudo-eigenvector \'^,E) includes at least one elliptic channel on at least one side. 
Therefore {n\'^,E) is not going to zero for n -> cx) or n — oo but {n\'i',E) is 
bounded. 



4. Normal forms of the transfer matrices 



For a formal eigenvector I^J^, E) of H, the coefficients are related by the transfer 
matrix and one has 7q^ (-ip*^) ) ~ ( (i-i^^g) ^ ■ Using the notations as in 



(3.111 and (3.121 one obtains from (3.11, (3.6) and (3.101 that 



'0,L 



E 



a<s 

.Ge{+,-} 



E 

.ee{+,-} 



(2sin(A:a))i/2 \ip^e 



E 



ee{+,-} 



•\/sinh(7) 



-e7o 



5© 



(2sin(/c„))i/2 y^^e 



-Qika 



E 



+ ,e 



^ sinh(7) 
ee{+,-} 



-070 



Working in the conjugate symplectic group one can diagonalize the hyperbolic chan- 
nels. In order to do this we define 



(4.1) C/=(^i,...,^Ar) e U(iV), 

(4.2) fc = diag(fci, . . . ,/cs), 7 = diag(7s+i, . . . , 7jv), m = diag(Ms+i, 
and the conjugate symplectic 2N x 2N matrix 

/ (sin(fe))"3 

(4.3) M 



.,Un), 



U 

u 





M(2sinh(7))"5 (2sinh(7))"5 

cos(fe)(sin(fe))^^ (sin(fc))5 

e"^(2sinh(7))"2 ue"' (2 sinh(7)) " 3 



\ e"^(2sinh(7)) 2 

Then M transforms the free transfer matrix to its symplectic normal form 



(4.4) 



'cos(fe) -sin(fc) 

Me"' 

sin(A:) cos{k) 

^ 
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and one obtains 

(4.5) ( 



a++a- \ 




( \ 


2a+ \ 






(a--a+) I 


= 


t(b--b+) I 


2a" / 




\ 2b- ' 



where T'' = Ar%%M . 



Note that T e Sp(2A^, C). To diagonahze the eUiptic channels for 7" ■ we need 



to conjugate M ^T^'^M by the Cayley matrix as defined in (2.10). This way we 



■■N 



e'" ~ 

tiCosh(7) usinh(7) 

e"''= 

sinh(7) tiCosh(7), 



obtain the normal form of the free transfer matrix in the Lorentz group U(A'^, N), 

Furthermore one obtains for CnT^C*^ G U(A^, N) that 
(4.6) 

5. Reduced transfer matrix 



N I (--AT 




We want to define a reduced transfer matrix relating the coefficients for the 
elliptic channels appearing in scattering states. This means we look for solutions of 
the equations above where oT = = 0. Given a^, and oT — the question is 
whether there exists a unique such that 6+ = 0. This is the case if the following 
(N — s) X (N — s) matrix 

/ 0, 



( 



■-(JV-ajxCN-s) "(JV-s)> 



-■(JV-s) X (N-s) 



(5.1) Ae 

is invertible. The indices indicate the size of the matrices. 



j-E 



5X(lV-a) 
,x(N-s) 



0, 



J(N-s)x(N-s) , 



Lemma 5.1. The matrix Ae is invertible for all but finitely many energies E in 
the spectrum of Hq . 

Proof. Let / be a bounded energy interval without parabolic channels. In I the 



elliptic and hyperbolic channels as well as the matrices U and u (as defined in (4.1 ) 
and (4.2|) stay the same. Now 7^^^ is of the form 

'E^ 




OX 



P{E) 



where P{E) is a polynomial in E of degree L — 1. Hence, we obtain from (4.31, 



( |4.5[ ) and (jSAj) for e / that 

(5.2) y/2 smh{j)e'"' Ae v/2sinh(7) = E^+Pi{E)e"'+e-'' P2iE)+e~'< P3iE)e''' 

where Pi{E), P2{E), P^iE) £ Mat(iV - s,C) are aU polynomials in E of degree 
L — 1. Letting A = diag(As+i, . . 

(5.3) 



, Xpf) one obtains from (3.9) that 



e-^ = e-^(^) = i 



i{A -E)- v/[w(A-£;)]2-4 



Since w(A - E) ^ 2cosh(7) > 2 for £; e /, the functions E ^ e^^t^) and E Ae 
can be extended to complex analytic functions on the strip / x iM c C. If the 



imaginary part '^(E) tends to oo, then e~'^ tends to zero. Multiplying (5.2 1 by 
E~^ and letting 3(i?) — >■ oo, the right hand side converges to 1. Therefore, Ae is 
invertible for large Q{E) and det(^£;) is analytic and not identical to zero. Hence, 
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det{AE) except for finitely many energies E in the precompact interval /. As 
there are only finitely many energies with parabolic channels, this shows the claim. 
□ 

Remark. An interesting question might be the meaning if Ae is not invertible. 
In this case Ae has a kernel and one can find a+ such that Aeo.^ — 0. This means 




23 + 







1 \ 






)- 


, i(b--b+) 1 




V 2b- / 



If one finds furthermore that 6+ = and 6^ = 0, then this corresponds to an 
eigenvector of H described by a+ and b~ and E is an eigenvalue of H . If the latter 
is not the case then we find a scattering state that has no elliptic channel on the 
left since a+ = = 0. From the interpretation of scattering states which will be 
given by the Lippmann Schwinger equation this means that there is an extended 
state or wave which is totally reflected. As we have seen, this happens only for 
finitely many energies. In particular, H has only finitely many eigenvalues. 



Let us now consider an energy E where is invertible. Then any vectors a' 



define a unique scattering state characterized by the coefficients a' 



,b+,b- 



letting 



one obtains 



b as defined in (3.11) and (3.12). More precisely, choosing vectors 



and 
and 



2a^ 



— a 






and has found all coefficients for the scattering state 
reduced 2s x 2s transfer matrix by 

(5.4) ^ 



In this case we define the 



a ' + a 

i(ar — a+) 

Another way to write would be 

(5.5) r^ = (Sg?g)r^ 



5+ 
i{b- 



-6^ 



^E 



1 



1 





The size of the first matrix on the right hand side of the equation is 2s x 2A^, the 
columns are divided into two blocks, each of size s, and the rows are divided in 4 
blocks of sizes s, N — s, s and iV — s in that order. This matrix is the same as the 



matrix O in Theorem 2.1 (iii). The last matrix is the transpose of the first one. 



2.1 



(iii). A conjugation of (5.4) by the Cayley matrix 



The matrix to the left and right of A^ is the same one that appears in (|5.1|) and it 
is equal to O* as in Theorem 
yields 

(5.6) ^"^ 



— i ' 



for r^ = Csr^ct 



Note, if s = then there is no hyperbolic channel and therefore all formal 
eigenvectors of H are scattering states and already relates the elliptic channels. 
Therefore, in this case one simply defines = ■ Then = CnT^C*^ and 



equations (4.5) and (5.4) as well as (4.6) and (5.6) are the same. In particular, T 
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is conjugate symplectic and is pseudo-unitary. This is actually always true, if 
the reduced transfer matrix exists. 

Proposition 5.2. The reduced transfer matrix is conjugate symplectic, i.e. G 
Sp(2s), and consequently, € U(s,s). 

Proof. Let Xi,yi G C, i = 1,2 and define Xi,yi by 



( J^M = I ^' ) , for i = l,2, then 3a^,6^ G C'^"" : T 



/ \ h. 



and one obtains 

As this is true for arbitrary Xi,yi one has {T^)* JsT^ — Js and hence is 
conjugate symplectic. □ 



= S:\y2 


- VlX2 








_ 1 ai 




- \ yi 




\ 



6. Scattering operator and scattering matrix 

Let n± = s— limt_>.;poo 6**^6"'*^" be the M0ller operators and S = 
the scattering operator. S commutes with Hq and can therefore be defined as 
operator on the energy shells for almost all energies E in the spectrum of Hq. 
So let |\E'i'„,i5) be some pseudo-eigenvector of Hq with such an energy E, then 
l^out^ ^) = '^l^hi' ^) is defined and also a pseudo-eigenvector of Hq with the same 
energy E. The subscripts 'in' and 'out' correspond to the physics intuition that the 
scattering operator maps the incoming states to the outgoing states. Furthermore, 
for almost all energies E the M0ller operators can be defined as maps from the 
energy shell with energy E with respect to the operator Hq , to the energy shell 
with the same energy E, with respect to the operator H . In this sense, we have 

(6.1) \^,E)=^l+\^l,E)^n-\^l^„E) , \^l^„E)^S\^l,E) , 

where E) is pseudo-eigenvector of H . One obtains from the Lippmann-Schwinger 
equations pUl [T5] . 

(6.2) {n\^,E)^{n\^l^,E) + Xim {n\{E - Hq - ie)-\H - Hq)\^ ,E) , 

e\,0 

(6.3) {7i\^>,E)^{n\^>l^„E) + ymi{n\(E-HQ-ie)-\H-HQ)\^,E) . 
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Inserting the partition of unity (3.3) and changing to an integral over the unit circle 



in the complex plane by substituting e*'^ = z, dk — i ^ dz gives 



{n\iE + Ho-te)-\H-Ho)\^,E) 

N 



J2 / dk{{n\^l,k){^lk\{E-H,-^e)-\H^Ho)\^,E)} 
„=i J--^ 



N 



= E 



dk 



ikn 



^ J_^ 27r 1 + le + e**^ + e-**^ - 



L-1 



J2e-'''^Vo.iVrn-W){m\^,E) 



m=0 



E 



dz 



E 



(fa Z- 



i_i 27ri 1 ^„ ^2 + + le - Aa) + 1 



¥':(Vm-M^)(m|vI/,i?) 



E E ..ewi - W){m\^, E), 



where is the solution of z^ + + — Aq) + 1 = which is inside the unit disc. 
Let a > s, then the solutions for e = are UaC^^" and ^„ converges to UQ,e~'*'° for 
e ^ 0. For the elliptic channels a < s the solutions for e = are e^**^" being both 
on the unit circle. As < ^ sign of its imaginary part is different to the sign 
of the imaginary part of + (i^^)"^ = — E — le. Since k^ & (0, tt) this leads 
to lime\__o — 6*'^° ^^'^ linie/'o £,a — 6^**^" . Hcncc by the calculations above and 
K2\, (lol) we get 



(6.4) 



= {n\^l.E) + E E Linfe) f*o.{V,r. ~ W){m\^ , E) 

a<s m—0 

+ E E .... - ^) 



a>s m—0 



(n|*,£;) = (n|<,„i?) + EE 



(6.5) 



EE 

a>s m—0 



2z sin(A;a) 



m—0 

e-7a _ e7» 



-ipliV„,-W){m\^',E) 



Thus, we see that in the hyperbolic channels, the extended state \'^,E) has only 
exponential decaying parts which justifies the definition for scattering states. More- 
over, ii\'^,E) is the scattering state associated to the coefhcients a+, a^, 6+, b~ and 
a+j &~ as in (3.11 ) and ( |3.12[ ), then for n — > zbcxi the equations give 



(6.6) 



\^l,E) = Y,[a+\^lE, 



6^1*° E - 



and 



■7) 



— ikn/ L 
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Therefore, the scattering operator reduced to the energy shell can be described by 
the 2s X 2s matrix defined by 



In particular, from (5.6) we obtain that T — CsT C* £ U(s, s) represents the S- 



transfer matrix. The existence of follows from (5.61 and the following theorem 
which is proved in Appendix [A} 

Theorem 6.1. For any matrix T € U(s, s) there is a unique unitary matrix S G 
U(2s) with the property that for any a+, a~, 6^ G one has 



(6.9) 



Putting equations (6.1), (6.6), (6.7| and (6.8) together one can write the op- 



erator iS as an integral over the energy E. The number of elliptic channels is a 
step function s{E). So far we considered one fixed energy and set the elliptic chan- 
nels to be the ones for a ~ 1, . . . , s. But when varying E one should take into 
account that the channels which are elliptic are different ones for different energy 
intervals. Therefore let a{E, 1), . . . , a{E, s{E)) denote the elliptic channels for E. 
Correspondingly for pseudo-eigenstates satisfying iS|^'['j^,i?) — |\E'°^t,£'), define the 



Then the scattering matrix satisfies (6.8 1 with 



coefficients a and b 

^ )^ and the analogue definitions for b^. Furthermore, 



a 



(a± 



■ '"a(B,s(E))J 



let eE.i,+ be the z-th and eE,i.- be the {s{E)+i)-t\i canonical basis vector of C^''*^'^^ 
for i = l,...,s{E). By (|6J]) and the matrix element e|j,,.^5^e£; 



corresponds to the contribution of a .^ to a 



The meaning of the other 
matrix elements can also be read off these equations and one finally obtains the 
following. 



Proposition 6.2. 

(6.10) 



The scattering operator S is given by 



S = dE 



J2 e^'^(^-^^-'--e^«)|vI/O(^,^.),i?,-0) {e%^^^^S^eE,^,e) {K(e,),E,Q\ 



i,j = l,....s(E) 
. e,0G{ + ,-} 

where the correction phase 9{E,i, j,Q,(Z)) is given by 

0{E,i,j,G,(Z)) = -SQ^^kE,iL - S0,^kEjL 

with 

E = -2cos{kEA) + K{E,i) , kEA e (0, tt). 



The phase 9{E,i,j,Q,0) comes from terms of the form 



appearing as factors 



in (6.6| and m.7l. Now we can finally prove Theorem 2.1 



Proof of Theorem 



2.1 



The direct integral / C^"*^^' dE is represented by 
functions f{E) with f{E) G C^^'^-* and the scalar product is given by {f\g) = 
J f{E)*g{E)dE. Let us define the unitary operator U : e(Z, C^) ^ /® C^'^^^dE 
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by 



s{E) 



.11) {U\n,j))iE) = E E 

*=i Ge{+,-} 



and the diagonal s{E) x s{E) matrix /cb by 

(6.12) ks = diag{kE,i, ■ . ■ ,kE,s{E)) 



Using (3.4) and (3.71 one obtains that 14 is unitary. The equations (6.10) and (3.8) 
can be written as 

(6.13) S =w\r (^^ ^'l'"^) 

(6.14) Ho I El■,_,|E^dE U . 



1 

e"**^^^ 



clE 



U 



A special case is H — Hq, where 5 = 1 and ~ (e'*E^ e'*"^^ ^ ^^^^ hence (6.13) 
gives U*U = 1. Equations (6.111 - (6.14) show Theorem 2.1 part (i). Part (ii) and 



(iii) follow from the equations (4.5), (5.1), (5.5), (5.6) and (6.^ 



7. AS LIMIT OF HIGHER DIMENSIONAL SCATTERING MATRICES 



In this section we prove Theorem 2.3 Recall that in the introduction we con- 



structed an ideal lead described by the operator Hj as in ( 2.17 1. The corresponding 
Hcrmitian matrix Wj E Her(Af) was defined by (2.18) which is equivalent to 



(7.1) 



Wi = C/diag(Ai,...,A„ £;,...,£;) (7* 



where U — (ipi, . . . , (^at) as in (4.1 ). The corresponding extended states of Hi as 
well as the wave-numbers ka for the energy E are the same as the once of Hq for 



a < s and given by (3.5) and (3.6). For a > s there are additional extended states 



of Hj for the energy E defined as in (3.6) with ka 



Inserting a piece of the cable Hq of length m followed by the scatterer and 
another piece of the cable of length m into the ideal lead Hj is described by the 
operator H{m) as defined in (2.19) (cf. Figure [T]) 



Similar to above one can introduce the vectors and in describing a 
formal solution (^„)„ of the eigenvalue equation H{m)'^ = E'^ for n < —m and 
n > L -t- TO, but this time there are only elliptic channels. 

The transfer matrix of the inserted piece given by the block of H{m) from n = 
-m to n = L + m - 1 is given by r{m) = {T°^^)"'TiL(X°^'^T , with T' 



as 



defined in (2.7). To get the relation between of and bj we have to follow the same 
steps as in Sections |4] and [S] Hence, let us introduce the matrix Mj similar to M 
in (|43| by 



(7.2) 



M, = 




(sin(fc))-i/2 O' 

10 

cos(fc)(sin(/£))"i/^ (sin(/c))i/^ 

1/ 
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tions as in (4.5), i.e. 



Let a^,a^ ,b^, satisfy the rela- 



a' +a 
2a+ 

2a" 




Then using (4.4) one obtains 
(7.3) M-^T{m)M 




i(e-""'=b--e*™'=()+~i 



Furthermore by (|4.3|) and the definition of Mj one has 
(7.4) 



10 

_ ( M(2sinh7)^2 (2sinh7) 



10 

e"^(2sinh7)~ 2 we^* (2 sinh 7) " 2 



(7.5) 



(7.6) 



To simplify notations let us define 

= 1 ( u 
V2sinh7 Ve""^ 

and 

J y/2 sinh7 



ue • 



ue 



Then the equations (7.3), (7.4), (7.5) and (7.6) yield 



(7.7) 



Cn MY^T{m)MiC 




e-'^'^b- 



By (2.15) the matrix on the left hand side of (7.7) is equal to the S'-transfer matrix 
describing the scattering of H{m) with respect to Hi. Let Sf{m) be the related 
scattering matrix and let us also introduce a phase normalization and consider the 
matrices 



Sf{m)^D"'Sf{m)D"', where D 



r"" 

M 

e"*'' 
m 



Then (7.7) and the relation between scattering and S'-transfer matrix (|2.8|) yield 
Sf{m) 



(7 



b+ 
b++ib- 



As the unitary group is compact, there is at least one limit point of this sequence, 
let us call such a limit point Sf. If 7 is a multiple of the unit matrix, i.e. 7 = 7!, 
then multiplying (7.8) by (2sinh7)5e~^™ and taking the limit m — >■ 00 along a 



sequence where Sf{m) converges to Sf{m) we obtain from (7.5) and (7.6) that 

(7.9) Sf 

As the reduced transfer matrix is supposed to exist, we can always choose to 
get any 6+ we want and hence, a~ and b'^ can be chosen independently. 



/ ^ 




( " \ 


/ (l+?ue^)d 




/ (1 — li/e"' )d \ 
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If 7 is not a multiple of the unit matrix, then for each 7^ we take the entries 
of oT and 5^ to be zero which belong to a 7,3 greater than 7^. Then we multiply 
(7.8) by (2sinh7Q,)3e~^°™ and take the limit m — 00. Doing this for any 7^ we 
also obtain ( |7.9| ) for any vectors a~ , 5+ by linearity. 

Furthermore one can choose a~ = and tune a+ such that 6+ = 0. Then and 
converge both to and a*, are related by the scattering matrix . Hence, 
in this case the limit m — )■ 00 of (7.8) along an appropriate subsequence yields 

1+^ 



(7.10) 



Sf 



E 











S 



E 



These two equations, (7.9) and (7.101, determine any limit point oiSf{m) uniquely. 



Therefore, the limit Sf = lim„ 
sf and given by 



, Sf (m) exists and there is a relation between 



(7.11) 

where 
(7.12) 



S' 



r 

R' 



Sf 



R 





T' 
















T' 





R' 

















[sinh(7)+zw] [cosh(7)]"^ G U(7V - 



This shows Theorem 12.3 



Appendix A. Proof of Theorem 16.1 



In this appendix we prove Theorem 6.1 It states that for any matrix T G U(s, s) 



there is a unitary matrix S £ U(2s) such that 

Let T — (p £). We first prove the existence and then uniqueness of S*. As 
T G U(s, s) one has 

A* A = 1 + C*C , AA* = 1 + BB* 

(A.2) D*D = 1 + B*B , DD* 1 + CC* 

A*B = C*D , AC* = BZ)* . 

As A* A > 1 there exists a unitary matrix and a real diagonal matrix Q > 1 
such that A* A = Ul^QUi^+. Define Ur,+ by 

(A.3) = AUl+^Q^ ^ A = Ur^+^Ui,+ 



Then U*^+Ur.+ = VQ^t^/,+^*-4C/;*+ 1 and hence Ur,+ is unitary. Fur- 
thermore one has C*C ~ A* A-1 = Ui ^{Q—l)Ui^+. Hence there exists Ur- G U(s) 
such that 



(A.4) Ur,-^Q-1 - ^ C = Ur.-^Q-lUi^+ 

Ur.- is uniquely determined if Q — 1 is invertible, otherwise it is not. Now define 

Ul- by 

(A.5) Ui^+ = ^/Q^U; _D ^ D = Ur,-^/QUi^- 
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Then using (A.2) and (A.4) one finds Ui-Ui_ ~ 1 and hence Ui- is also unitary. 
Furthermore one obtains using (A.2), (A.3), (A.4) and (A. 5 1 

B ^ {A*)-'^C*D 



(A.6) 



Now using (A.3), (A.4), (A.5) and (A.6 1 one obtains 

\{ VQ VQ^ 

Ur,-) \VQ^ VQ 



(A.7) 



r 



Ur,+ 





Ui.- 



As Q > 1, one has 1 — Q ^ > and hence -^/l — is a weU-defincd, non- negative 
diagonal matrix. Thus we can define the unitary Matrix 

-Vl-Q-' VQ^ \ fUi.+ 

u: 



(A 



5* 









Ur.,+ 



It is now easy to check, that S fulfills ( A.l ) 



To prove uniqueness of S assume S also fulfills ( A.l ). But then ( A.l ) implies for 
any vector v E C^'* that Sv = Sv and hence S ^ S. 



□ 



Remark. The converse is not true. One cannot find a matrix T for all unitary 
matrices S G U(2s) such that the relation above is fulfilled. Looking at the block 
structure S = ( ^ ) The matrix T exists if T is invertible (which is equivalent to 
T' being invertible). T and T' are related to the transfer of waves. If they are not 
invertible, then there is one planar wave which is totally reflected by the scatterer. 
Hence a transfer does not occur for this wave and the transfer matrix is not defined. 
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